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Abstract 

We present an investigation of the atomic and electronic structure of graphene monolayer islands 
on the 6iif-SiC(0001)(3 x 3) (SiC(3 x 3)) surface reconstruction using scanning tunneling microscopy 
(STM) and spectroscopy (STS). The orientation of the graphene lattice changes from one island 
to the other. In the STM images, this rotational disorder gives rise to various superlattices with 
periods in the nm range. We show that those superlattices are moire patterns (MPs) and we 
correlate their apparent height with the stacking at the graphene/SiC(3 x 3) interface. The contrast 
of the MP in STM images corresponds to a small topographic modulation (by typically 0.2 A) of 
the graphene layer. Prom STS measurements we find that the substrate surface presents a 1.5 
eV wide bandgap encompassing the Permi level. This substrate surface bandgap subsists below 
the graphene plane. The tunneling spectra are spatially homogeneous on the islands within the 
substrate surface gap, which shows that the MPs do not impact the low energy electronic structure 
of graphene. We conclude that the SiC(3 x 3) reconstruction efficiently passivates the substrate 
surface and that the properties of the graphene layer which grows on top of it should be similar to 
those of the ideal material. 



I. INTRODUCTION 



Fascinating properties have been predicted and observed for monolayer graphend^^. 
Among them one finds the anomalous quantum Hall effeclP^, the Klein tunneling 
phenomenorP^, and weak (anti) localization effects^'^. Moreover, suspended graphene shows 
exceptionally high carriers 

mobility™ 

even near room temperatur^^. These features make 
graphene an attractive material for the investigation of original physical phenomena?2II3! and 
for the development of devices such as transistor d^^*^^ ! and captor^. 

The physical properties of free standing graphene are intimately linked to the presence 
of two equivalent carbon sublattices commonly called A and B. Usually, graphene layers are 
supported on a substrate and the interaction between the electronic states of the substrate 
surface and the pz orbitals of the C atoms can significantly alter the electronic structure -and 
thus the properties- of the material. This has been shown recently by angle resolved pho- 
toemission for graphene elaborated on metal surfaces where this coupling modifies the band 
dispersion close to the Dirac poinlP^^*^^*^, suppressing the "Dirac cones" . The investigation 
of the atomic and electronic structure of the interface between graphene and the substrate is 
thus of primary importance. This is in particular the case for few layers graphene grown on 
SiC substrates, where as-grown samples are used for physical measurement since the 

doped graphene layers close to the interface should give the largest contribution to electrical 
transport!^. 

Few layers graphene are obtained by high temperature treatment of the polar faces of 
SiC substrate j^^ * ^^ * ^^ ^ Usually commercial hexagonal {AH or QH) substrates are used. They 
have two different faces, the (0001) one (the Si face) and the (0001) faces (the C face). The 
interface between the Si face and the graphene overlayer has been extensively studied in the 
last few years. The current model for this interface is that the first graphitic layer strongly 
interacts with the substrate, giving rise to the x Qy/3)R30° {QR3) reconstructiorP^'^. 

Covalent bonds form between Si atoms of the substrate surface and the graphene layer, which 
results in the suppression of the Dirac cones characteristic of graphen^^^'^^'^. This model is 
supported by photoemission datsP^. Accordingly no graphene contrast has been detected in 
STM images of the 6R3 reconstructiorJ^ ^^ * ^^ * ^^ !, which is usually called the "buffer layer". 
The electronic structure of graphene is developed only for the second C plan^^MZMES]^ where 
a band structure very similar to the Dirac cones has been observed experimentallj^^^'^. The 



2 



question of a possible perturbation of the electronic structure of the graphene layer due to 
an interaction with the buffer layer remains operP^^^. Nevertheless, the honeycomb contrast 
expected for ideal graphene is observed by STM on this second C pland^^ ^ * ^^ * ^^ l Moreover 
the analysis of the standing wave patterns indicates that the electronic chirality of graphene 
is preservecP^. 

The interface between graphene and the C face has been less extensively studied. It has 
long been known that the growth is quite different on the C and the Si face^-. Graphitic films 
grown in UHV conditions on the C face exhibit some rotational disordei'^^'^. This disorder 
already exists for the first C layei'^^'^^^. Interestingly, it was found using photoemission that 
the interaction between the first C layer and the substrate was much weaker than for the Si 
face: no buffer layer is detected in core level spectroscop ^^" ' ^^' and the band structure of this 
layeil^Sl resembles the one of graphene. The situation is however complicated by the facts that 
i) two different pristine reconstructions of the substrate -the SiC(2 x 2)^ and the SiC(3 x 3)- 
exist at the interface below the graphene layei'^^'^^^ and ii) that several orientations exist 
for the graphene islands for each reconstruction, leading to different superlatticed^. A 
systematic analysis of the interface for the two different substrate reconstructions aiming 
at understanding their atomic and electronic structure for the different orientations of the 
graphene layer is thus needed. This is best achieved by STM, which can address each 
individual island. 

In a previous paper we have shown that a graphitic signal could be observed at low bias 
for both the SiC(2 x 2)c and the SiC(3 x 3) reconstructions, indicating a weaker coupling 
with the substrate than on the Si fac^^. A recent ab-initio calculation has shown that the 
reduced interaction in the case of the SiC(2 x 2)c reconstruction is due to a passivation of 
the substrate surface by Si adatoms^^. The linear dispersion of the graphene bands close 
to the Dirac point is preserved, but a residual coupling with the substrate, also evidenced 
by STM, was found. In the present paper we concentrate on the graphene islands formed 
on the SiC(3 x 3) interface reconstruction, which are called GS x 3 islands afterward, 
where the interaction with the substrate seems to be even smallei'^. We first analyze the 
geometric structure of the superlattices. We show that they are moire patterns and we 
relate their apparent height to the local stacking at the interface. We then analyze the 
electronic structure of the GJi x 3 islands compared to that of the bare substrate SiC(3 x 3) 
reconstruction. A wide surface bandgap (of width 1.5 eV) is found by STS in the electronic 
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structure of the bare reconstruction, which persists below the graphene layer. The Fermi 
level of graphene is located in the vicinity of the top of this surface gap. In STM images 
a graphene signal dominates inside the substrate surface gap, and the STS data explain 
the high bias "transparency" of graphene. A comparison between STM images of G-3 x 3 
islands for a specific orientation with previous ab-initio calculations -as well as with the 
case of the Si face- indicates that the substrate reconstruction is responsible for the weak 
graphene-substrate interaction. Finally we show that the moire pattern is essentially of 
topographic origin. It is associated with small undulations of the graphene layer. From STS, 
these undulations do not lead to heterogeneities in the electronic structure of graphene, at 
variance with the case of more strongly interacting systems such as graphene on Ru(OOOl)'^. 
From these data we conclude that the GJi x 3 islands should be a system close to ideal, 
uncoupled, graphene. At present it is not clear whether the GJi x 3 structure is present 
at the interface for few layers graphene films elaborated at high temperature in non UHV 
conditiond^. Our results indicate anyway that manipulating the atomic structure of the 
surface can be a useful way to modify the coupling at the interface, as shown previously for 
metal substrate^. 

II. EXPERIMENT 

The sample preparation and characterization were conducted under ultrahigh vacuum. 
The sample graphitization was performed in-situ by following the procedure presented 
previously"*-. The n doped 6if-SiC(0001) samples were first cleaned by annealing under 
a Si flux at 850°C. The already reported SiC(3 x 3) reconstruction'*'^ was obtained by further 
annealing at 950 — 1000°C. After annealing steps at increasing temperature, a graphene cov- 
erage of less than a monolayer is finally detected by low energy electron diffraction (LEED) 
and Auger spectroscopy. At this stage, the LEED patterns show SiC(3 x 3) and SiC(2 x 2) 
spots and a ring-shaped graphitic signal with modulated intensity* ^° * ** * *^ l 

The STM and STS measurements were made at room temperature with mechanically 
cut Ptir tips. 5 samples were investigated, using more than 10 macroscopically different 
tips. The samples morphology observed by STM was similar to previous results, with the 
presence of bare SiC(3 x 3) reconstructed substrate domains, graphene monolayer islands on 
the SiC(3 x 3) reconstruction (G_3 x 3) and on the SiC(2 x 2)c reconstruction (G_2 x 2) and 
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also few multilayer island^^^^^. The focus of this paper is the structure of G_3 x 3 islands 
and several dozens of them were observed. 

III. RESULTS AND DISCUSSION 

A. Superlattices and local stacking of monolayer graphene on the SiC(3 x 3) re- 
construction 

STM images show that GJi x 3 islands present superlattices (SLs) of various periodicities 
in the nanometer range (see Fig. [T] (a)). Contrary to graphene on SiC(OOOl) (Si face), 
graphene on SiC(OOOl) (C face) exhibits a significant rotational disorder, already from the 
first graphene layer. This results in a ring-shaped graphitic signal on LEED patterns (Fig. 
[l] (b)). A previous STM study established that the SL period depends on the orientation 
angle of the graphene layer with respect to the substrate surface lattice^^. In this section 
we present a more quantitative analysis of the SLs that identifies them as moire patterns 
(MPs). 

Moire patterns arise from a non linear composition of two periodic latticed^. They appear 
as an additional periodic lattice of larger period than the two components. For example, 
MPs are observed by STM on graphite^^ and few layer graphene samples with rotational 
stacking fault^. They also show up when two lattices of different lattice parameters are 
superimposed as for graphene monolayer on transition metal^^ ^^ * ^^ l For two periodic lattices 
with reciprocal lattice vector ki and k2 respectively, the resulting MP is characterized by 
the reciprocal lattice vectoil^SEB 

kM = ka - ki. (1) 

In the system we consider, the SiC(3 x 3) lattice parameter being almost 4 times bigger 
than the one of graphene, high order spectral components have to be considered. As we 
can see on the LEED pattern in Fig. [T] (b), first order SiC(3 x 3) spots are located far 
away from the graphitic signal. According to equation (1), moire patterns constructed on 
these spots and any graphene spot would have a smaller period than the SiC(3 x 3) lattice, 
which cannot explain the observed SLs. SiC(l x 1) spots can also be ruled out for similar 
reasons. From the LEED pattern of fig. [l]b), the reciprocal lattice vectors of the SiC(3 x 3) 
reconstruction most likely to lead to MPs with periods in the nanometer range are the high 
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order (4, 0), (3, 1), (2, 2) ones and their symmetric counterparts in the reciprocal space. Fig. 
[T] (c) provides an illustration of a MP construction associated to the (2, 2) SiC(3 x 3) spot 
of the LEED pattern, for a graphene island of orientation a with respect to the SiC surface 
lattice. 

Thus, we calculated the moire periodicity P as a function of the graphene orientation 
angle a with respect to the SiC surface lattice with a ranging from 0° to 30° (due to the 
symmetry of the system seen on the LEED pattern). For each of the three relevant SiC(3 x 3) 
Fourier components, we use equation Q and P{a) = 2it /{kM (a) cos{n / 6)). The three 
resulting P{a) curves are plotted in Fig. [l] (d). We have also measured moire periodicities 
versus graphene orientation angles on STM images of monolayer GS x 3 islands (such as 
Fig. [l] (a)), with an accuracy of ± 0.1 nm and ± 1 ° respectively. As represented on Fig. 
[T] (d) , experimental data do fit very well with calculations. For a given angle, the largest 
period - which corresponds to the best match in reciprocal space - is generally predominant 
in the images. We also note that most studied islands exhibit an orientation angle between 
15° and 30°. This is consistent with the peculiar rotation angle distribution revealed by 
LEED'^^J^. We will thus concentrate on these orientations in the following. To summarize, 
we interpret superlattices on G_3 x 3 monolayer islands as high order MPs, resulting from 
the superposition of the SiC (3 x 3) and the graphene-like lattices. Note however that the 
moire interpretation is essentially geometric and that it does not give any information on 
the nature of the interaction between graphene and its substrate. We shall consider this 
point in section [III C 



We now focus on the atomic structure and on the stacking for graphene islands with a 
MP constructed on the (2, 2) and (3, 1) SiC(3 x 3) Fourier components which are the most 
common on our samples (15° < a < 30°). As shown on Fig. [l] (d), the corresponding moire 
periodicity is maximum for a graphene orientation angle a of 30° and 13, 9° respectively. 
Low bias STM images (see Fig. [2] (a), (b)) show that the MPs observed for angles close 
to these two values exhibit inverted contrasts: "ball-like" for a close to 30°, "hole-like" for 
a close to 14°. At the atomic scale, a well defined honeycomb pattern characteristic of 
monolayer graphene is observed at low bias for both orientations (see Fig. |2](a), (b)). We 
point out that the moire contrast on G'_3 x 3 islands shows no variations with the tip and 
tunneling conditions: bright areas remain bright at any bias and for all tips tested (see Fig. 
[5l). 
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In order to understand the variation of the MP contrast with angle a, we studied the 
local stacking of the graphene and SiC(3 x 3) lattices for a close to 30° and 14°. As already 
mentioned in previous paper^^^^^^ISIli]^ graphene appears transparent on high bias STM im- 
ages so that the interface -the SiC(3 x 3) reconstruction in the present case- becomes visible. 
Conversely, atomic resolution on graphene is obtained on low bias images. Thus, stacking 
can be observed using two different approaches: by dual bias imaging at low and high bias or 
by imaging at an intermediate tunnel bias voltage, which corresponds to a crossover between 



these two extreme situations (to be discussed in section IllB). The latter type of image is 
represented in Fig. ^c) and [2](d) for islands with a = 29° and 15° respectively. For this 
sample bias {Vs = —1.65 V), the graphene and SiC(3 x 3) lattices appear simultaneously. 
Schematic reproductions of the images in Fig. |2]^e) and (f) give a clear view of the local 
stacking (the same result is found from dual bias imaging). Since no established structure 
model for the SiC(3 x 3) reconstruction exist, we represent in Fig. [2] (e), (f) the states 
detected by STM on the substrate surfac^^ filled (empty) states are represented in red 
(light-red) (dark gray and light gray in the printed version). 

In the a = 30° case (Fig [2] (e)), the graphene and SiC(3 x 3) lattices are quasi com- 
mensurate. The corresponding common Wigner-Seitz cell is represented by solid lines. The 
apparent height of the MP is maximum in the center of the cell and minimum on its edges. 
These areas correspond to two different types of stacking. At the center of the cell (circled 
area), SiC(3 x 3) states are located under the center of graphene hexagons (i.e. no C atom 
is in coincidence with SiC(3 x 3) states). On the edges (dark area), all SiC(3 x 3) states 
have C atoms or C-C bonds on top of them. 

In the a = 15° case, the moire corrugation is inverted. The apparent height of the moire 
is minimal at the center of the cell and maximal at its edges. Now the stacking at the center 
of the cell (circled area) is such that every SiC(3 x 3) state has C atoms or C-C bonds 
directly above which is similar to the stacking in the dark area of the a ~ 30° case. At 
the edges of the cell, a significant amount of SiC(3 x 3) states are located under the center 
of graphene hexagons, as for the bright regions in the a ~ 30° case. Therefore, the local 
stacking of bright (high) and dark (low) areas is the same for the two kinds of MP contrast, 
"ball-like" {a ^ 30°) and " hole-like" (a ^ 14°). The apparent MP "contrast inversion" 
arises from changes in the local stacking induced by the graphene rotation. 

Islands with a ~ 30° deserve particular attention because they allow a direct comparison 
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with the experimental results for the Si face and with theoretical calculations. For a = 30°, 
the graphene and the SiC(l x 1) lattices are (quasi)-commensurate with a 6^/3x 6\/3R{30°)- 
SiC (6i?3) common cell (or a (13 x 13) graphene cell). This is the configuration which is 
observed for the Si fac d^^ * ^^ * ^ , the layer orientation is imposed by the substrate and is 
therefore the same on the whole sample. A strong interaction between the first graphitic 
layer ( "buffer layer" ) and the substrate occurs so that only the second layer shows graphene 
propertie d^^ * ^^ * ^^ * ^ . In particular, no honeycomb contrast characteristic of graphene has ever 
been observed in STM studies of the 6R3 phase of the Si face -corresponding to the first 
graphitic layer or "buffer layer"- since it lacks vr states in the vicinity of the Fermi leveP^. 
Additionally, this 6R3 usually gives rise to a dominant SiC(6 x 6) superstructure in STM 
imaged^, although high resolution images reveal the actual 6-R3 periodicitj^. 

A totally different situation occurs for graphene on the SiC(3 x 3) reconstruction of the 
C face. For a = 30°, Fig. [3] (a), we actually observe a 2^/3 x 2-\/3_R(30°) with respect to the 
SiC(3 X 3), which corresponds to the actual 6R3 (and not SiC(6 x 6)) superstructure with 
respect to the SiC (1x1). More important, the honeycomb contrast of graphene clearly 
shows up at low bias (see inset in Fig. [3]). This demonstrates that the graphene states are 
present close to the Fermi level and thus implies a comparatively much weaker interaction 
with the substrate compared to the Si face. This weak coupling probably results from the 
presence of the SiC(3 x 3) surface reconstruction below the graphene layer. Indeed, ab-initio 
calculations performed for a graphitic C layer on the ideal (non-reconstructed) C face for 
this orientation {a = 30°) indicate a strong bonding to the substrat d^^ * ^^ * ^^ ! and subsequently 
the disappearance of the vr states at low energy, as for the Si face. This suggests that the 
SiC(3x3) reconstruction efficiently passivates the substrate surface for the C face, preventing 
the formation of chemical bonds with the graphene layer. Similar results were obtained from 
ab-initio calculations for graphene monolayer on the SiC (2 x 2)c reconstructiori^, a system 
that coexists with G_3 x 3 on the SiC(OOOl) graphitized surface. 

Another observation we have made on G-3 x 3 islands for a ~ 30° is that the orientation 
of the graphene layer is not locked to 30°. To see that, we took advantage of the fact 
that the MP orientation with respect to the SiC lattice is highly sensitive to the rotation 
angle a. Indeed, for a change by 1° of a, the orientation of the MP changes by 13° (for 
a ~ 30°). Thus, we could identify some slight deviations {Aa < 1°) from the quasi- 
commensurate {a = 30°) configuration. Fig. [31(b) provides an illustration of this effect: the 
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MP significantly differs from the 6R3 supercell although a remains close to 30° (we measure 
a = 29 ± 1°). The fact that such islands exist suggests that the 6R3 (quasi-commensurate) 
configuration does not lead to a notable energy reduction, at variance with the case of the 
Si face. This is consistent with the absence of covalent, directional bonds between graphene 
and substrate at the interface for GJi x 3 islands. 

B. Electronic structure of monolayer graphene on the SiC(3 x 3) reconstruction 

After these mainly structural considerations, we focus on the electronic structure of the 
graphene overlayer and of the SiC(3 x 3) reconstruction. We first compare the electronic 
structure of the G.3 x 3 and SiC(3 x 3) phases using current imaging tunneling spectroscopy 
(CITS). This technique consists in acquiring a constant current image and an I(V) curve 
after each of its points. For each spectrum, the feedback loop is turned off and the sample 
voltage (Vs) is ramped between preset values. I(V) curves are then numerically differentiated 
to get dl/dV conductance curves which are -in first approximation- proportional to the local 
density of states (LDOS) of the sample surface. In Fig. |4| we present CITS data acquired 
on a region (see insert in Fig. |4](a)) with the bare SiC(3 x 3) reconstructed substrate surface 
(right) and a G_3 x 3 island (left), so that both region are probed with the same tip. Fig. |4] 

(a) shows three I(V) curves, one for each type of surface, spatially averaged over the boxed 
regions (300 points each), and one for the edge of the graphene island (averaged over 15 
points). The dl/dV spectra for the x 3 island and the bare substrate are given in Fig. 

a(b)- 

For the bare SiC(3 x 3) reconstruction, the I(V) curve in Fig. |4] (a) displays a dramatic 
reduction of the current between Vs = —1.4V and +0.1V. This feature is still visible - 
although less marked- in the I(V) curve for the graphene island. The curve obtains on the 
edge of the island indicates that the lack of current at low bias does not arise from the 
electronic structure of the tip. These observations suggest the presence of a surface bandgap 
associated to the SiC(3 x 3) reconstruction that subsists under the graphene layer. However, 
a residual current related to in-gap states is detected in the surface gap of the bare SiC(3 x 3) 
reconstruction (between = —1.4V and —0.5V). To further study the electronic structure 
of the G-3 X 3 island and of the bare substrate, the conductance curves presented in Fig. |4] 

(b) are analysed in the following. 
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The SiC(3 x 3) spectrum exhibits a region of minimum conductance ranging from = 
—1.4 V to +0.1 V. These values are only weakly dependent (within 0.2 V) of the tip and 
sample. Hence, the SiC(3 x 3) reconstruction presents an asymmetric surface bandgap, 
with the Fermi level close to the bottom of the conduction band, as expected for a n- 
type semiconductor. A broad structure centered around Vs = —0.7 V is also detected. It 
is ascribed to in-gap states. An additional CITS study of the SiC(3 x 3) reconstruction 
suggests that they arise from a subsurface atomic layer (not shown). For the occupied 
states, these observations are consistent with Angle Resolved Photoemission Spectroscopy 
(ARPES) dataP^of the bare and lightly graphitized SiC(0001)(3 x 3) surface where a large 
intensity for binding energies larger than 1.5 eV and a residual emission between 0.5 eV and 
1.0 eV binding energies are detected. 

The GJi X 3 spectrum has similarities with the SiC(3 x 3) spectrum. In particular, the 
structure at Vs = +0.7 V (above the top of the substrate surface gap) and the rapid increase 
of conductance below the bottom of the substrate surface gap {Vs ~ —1.4 V) are clearly 
observed. These structures arise from the underlying SiC(3 x 3) reconstruction. The Fermi 
level of graphene (Ep in Figji] (c)) is located close to the top of the substrate surface gap. 
Inside the SiC(3 x 3) surface bandgap, an additional - though rather small - density of states 
originating from the graphene layer is detected. In other words, outside the (bare) SiC(3 x 3) 
surface bandgap, the signal is dominated by the contribution of the substrate which explains 
the transparency of graphene at high bia d^^ * ^^ * ^ ^*^. 

Note that the surface bandgap of the substrate remains unchanged below the graphene 
layer. This again suggests a weak graphene - substrate interaction since strong coupling 
would also affect the electronic states of the SiC(3 x 3) reconstruction. Moreover, for an 
extended energy range within the surface bandgap (from —1.4 eV to +0.1 eV) the density 
of interface states susceptible to interact with graphene states is quite small (from Fig. |4]). 
This is consistent with graphene-like atomic contrast on low bias STM images presented 
here (Fig. |2]) and in previous paper^. Nevertheless, moire patterns are still visible on the 
graphene islands at energies within the SiC(3 x 3) surface bandgap. This gives evidence for 
a residual effect of the substrate. Following experiments aim at discriminating between a 
topographic or an electronic effect for the MP contrast. 

In Fig. |5] (a) we present a series of STM images of the same area of a G'_3 x 3 island with 
a = 16° at various sample bias voltages. At Vs = —2.5 V, the SiC(3 x 3) reconstruction 
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is clearly visible while no evident graphene signal is detected. At Vs = —1.5 V, atomic 



(see Vs = —1.0 V and Vs = +10 mV panels). From Vs = +0.5 V to higher biases, no 
more evident graphene signal is visible and the SiC(3 x 3) signal reappears. Thus, atomic 
resolution on graphene is mostly achieved within the SiC(3 x 3) surface bandgap, as expected 
since the DOS arising from the SiC(3 x 3) is small within the surface bandgap. If we now 
concentrate on the moire pattern - the signal with period P = 4.5 nm in the images of Fig. 
[5] (a) - we notice that its corrugation is obviously much smaller at high biases than at low 
biases. 

To complete these observations, we have measured the peak-to-peak moire corrugation 
amplitude as a function of the sample bias voltage (see graph in Fig. |5](c)) using two different 
methods: i) profiles laterally averaged over 0.6nm on raw images or ii) profiles taken on low- 
pass filtered images in order to get rid of the high frequency atomic corrugation. Both 
methods gave the same results. The former method is illustrated in Fig. [s] (b) for a profile 
taken on the low bias image {Vs = +10 mV) of Fig. |5] (a). We stress that the measurements 
reported in Fig. [5] (c) were actually made at several spots on larger images. The uncertainty 
in the measurement is estimated to be of ±0.025 A . It results from the residual contribution 
of the atomic corrugation (SiC(3x3) at high bias or graphene at low bias as in Fig. |5](b)) and 
to some inhomogeneities in the corrugation of the moire patterns. The data presented in Fig. 
[5] (b) were acquired with different tips on two different GS x 3 islands (labeled "Island 1" 
and "Island 2") of approximately the same orientation: a = 16°, P = 4.5 nm (STM images 
in Fig. [5] (a) were obtained on "Island 1"). A typical SiC(3 x 3) average spectrum is also 
given on the graph in order to locate the substrate surface bandgap. The behavior is similar 
for both set of measurements with the following characteristics: within the SiC(3 x 3) surface 
bandgap, the moire corrugation amplitude is constant and equal to 0.25 ± 0.025 A. Outside 
the surface bandgap, when the SiC(3 x 3) contribution to the tunneling current becomes 
dominant, the corrugation dramatically decreases -or even vanishes. Other G'_3 x 3 islands 
with different orientations (including angles close to 30°) show similar behavior of the moire 
corrugation as a function of bias. These observations imply that the moire corrugation is 
associated to the graphene layer and not to the SiC(3 x 3) reconstructioEp^l. Moreover, since 
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the corrugation is independent of the bias in a large voltage range spanning the SiC(3 x 3) 
surface bandgap, it is most probably of topographic origin. The corrugation amplitude 
varies with the moire period, between 0.15 A and 0.25 A (smaller corrugations correspond 
to moire patterns of smaller periods). 

From the above observations, the graphene-substrate distance changes by 0.2 A from 
the highest to the lowest areas of the moire pattern. To look for a possible change in the 
electronic properties correlated to these soft "ripples" , we have performed scanning tunneling 
spectroscopy on several islands. In Fig. [6] we present CITS results on a GJi x 3 island with 
a ~ 30° (see insert) - similar results were obtained for other values of a between 15° an 
30°. The setpoint was chosen within the substrate surface bandgap {Vs = —1.0 V) with a 
setpoint current of 1.0 nA in order to probe essentially the graphene states. The graph shows 
two spectra, one corresponding to the highest regions of the moire pattern and the other 
to lowest regions. Both are averaged over 24 points. The spectra coincide from Vs = —1.0 
V to +0.15 V. This demonstrates that the electronic structure of the graphene overlayer is 
homogeneous and thus not affected by the moire pattern in a wide energy range spanning 
the Fermi level (and located in the substrate surface bandgap). Note however that the 
spectroscopy measurements are conducted at room temperature and the energy resolution 
is thus limited to ~ 0.1 eV--. From = +0.2 V to +0.5 V, more signal is detected on 
the low regions than on the high ones. This discrepancy arises from topographic effects 
(decrease in MP corrugation) discussed in connection with Fig. |5} 

Another important question is the position of the Dirac point. Previous ARPE^^ and 
transport measurement^ assess that it is located 0.2 eV below the Fermi level. But these 
techniques are non-local and give therefore an average value of the doping of the graphene 
layer. Importantly, the underlying reconstruction of the substrate is not identified in the 
probed region. STS is thus in principle the most adapted technique for answering this 
question. In Fig. |6} we find a rather structureless spectrum with a "flat" minimum ranging 
from Vs = —0.2 V to V. Some other spectra showed a well defined minimum located 
around Vs = —0.25 V. However, due to a significant variability in our measurements of the 
dl/dV curves between Vs = —0.5 V and V, we refrain from giving a definite value for the 
position of the Dirac point. 
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C. Graphene on 6i7-SiC(0001)(3 x 3): an almost ideal graphene layer? 



From measurements presented in Fig. [s} we find a topographic corrugation of 0.15 A 
to 0.25 A of the graphene monolayer while no long range topographic modulations were 
observed on the bare SiC(3 x 3) reconstruction. The period of the graphene topographic 
modulation is related to its orientation with respect to the substrate reconstruction and 
follows a moire model (discussed in connection with Figjl]). This means that the graphene 
corrugation is induced by the SiC(3 x 3) reconstruction. More precisely, the graphene - 
substrate distance is governed by the local stacking of the SiC(3 x 3) and the graphene 



overlayer as shown in section III A 



Such an effect has already been reported for graphene on transition metals. It is instruc- 
tive to compare our data with a well documented case of a relatively strong coupling, such 
as graphene on Ru(OOOl), where the vr bands of graphene are strongly perturbed by interac- 
tion with the substrat d^^ * ^^ l This system presents a moire pattern (P = 2.9 nm) caused by 
the lattice mismatch between graphene and Ru. A signal with the periodicity of graphene 
is observed by STM!^^*^ but the contrast changes from honeycomb in the high region to 
triangular in the low axead^^^^^. This is at variance with the uniform honeycomb pattern 
we observe on GJi x 3. DFT calculation^, surface X-ray diffraction^^ STMSMH and core 
level spectroscopjJ^ conclude that lower areas of the graphene layer strongly bond to the 
substrate. 

STS results presented for graphene on Ru(OOOl) in Ref. EH are of particular relevance 
for the purpose of our study. It shows dl/dV spectra with significant spatial variations 
correlated to the moire pattern. This finding was interpreted in Ref. HH using a generic 
model where a periodic potential -with the periodicity of the MP- is applied to a flat and 
isolated graphene layer. This leads to the LDOS modulations observed by STM, which also 
corresponds to charge inhomogeneities in the graphene layer. In Ref. EH it was inferred 
from DFT calculations that the spatial variations of the STS spectra should be attributed 
to the spatially heterogeneous bonding between graphene and Ru. This is clearly different 
from the behavior we observe for the dl/dV spectra on G'_3 x 3 (Fig. [6]), and we thus 
conclude that neither charge modulation nor local (periodic) bonds formation occurs in this 
system, whatever the orientation angle a. Incidentally, even for graphene on SiC(OOOl), 
where the graphene overlayer (second graphitic plane) is known to be well decoupled from 
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the substrate, spatial variations of the dl/dV spectra have been reported close to the Dirac 
poinlP^. Therefore graphene on the GJi x 3 islands may be quite close to ideal graphene, 
due to a weak interaction with the substrate reconstruction. 

We now briefly discuss the origin and the influence of the corrugation of the graphene 
layer which give rise to the MP. Since strong periodic bonding to the substrate can be ruled 
out from our data, these topographic modulations probably come from a weak, possibly Van 
der Waals-like, interaction that depends on the local stacking. Note that the corrugation 
we measure is small, typically 0.2 A Peak to Peak (PP) for wavelengths P in the range 
2 — 5 nm. The consequence of such "ripples" on the electronic structure of isolated graphene 
layers has been estimated in previous papers. For a graphene layer with a modulation of 
pseudo-period P' = 1.9 nm and an amplitude of 0.4 A PP, DFT calculationd^^*^ show no 
significant modification of the electronic properties with respect to the fiat configuration. 
In particular, it does not open a gap at the Dirac pointP^. Even on an isolated strongly 
corrugated monolayer (1.5 A PP for a period P ^ 3 nm), other ab-initio calculations shows 
that the LDOS of graphene remains linear within ±1 eV from the Dirac point in the high 
(and low) region^. Thus we believe that the small topographic corrugation we observe 
should have only a limited effect on the electronic structure of graphene close to the Dirac 
point for x 3 islands. However, experiments with an improved resolution should be 
performed to search for -or to rule out- a possible influence of the superperiod (MP) on the 
band structure of graphen d^^ * ^^ l 

IV. CONCLUSION 

We have investigated by STM and STS graphene monolayer islands grown under UHV on 
the SiC(3 x 3) reconstruction of the 6iJ-SiC(0001) surface. These islands present different 
orientations with respect to the substrate. From STM topographic images with atomic res- 
olution, we find that the various superstructures with periods in the nm range observed on 
the islands can be interpreted as moire patterns arising from the composition of graphene 
and high order SiC(3 x 3) lattice Fourier components. We show that the moire contrast 
corresponds to topographic modulations in the graphene layer of typically 0.2 A. The lo- 
cal graphene-substrate stacking in the low and high regions of the moire pattern could be 
obtained using the transparency of the graphene in high-bias images. Our STS study of 
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the SiC(3 x 3) substrate reconstruction reveals a surface bandgap (typically ranging from 
-1.4 eV to +0.1 eV) that persists under the graphene monolayer. This characteristic ex- 
plains the variations with sample bias voltage of the graphene/substrate signal ratio for 
STM topographic images and for STS. Further STS measurements show that the electronic 
structure is spatially homogeneous for any orientation of the graphene layer which indicates 
a weak graphene substrate interaction. This is confirmed by the absence of preferential 
graphene orientations even for an almost commensurate configuration {a ~ 30°). The weak 
interaction is achieved thanks to the surface reconstruction that efficiently passivates the 
SiC substrate since ab initio calculations for a bulk-truncated SiC(OOOl) surface have pre- 
dicted strong interaction and covalent bonds formation between the graphene layer and the 
substrat^2MHESl_ This suggests that the graphene-substrate coupling can be tuned using 
post-treatments that alter the substrate surface reconstruction. Finally, in an energy range 
of ±100 meV spanning the Fermi energy, very few substrate interface states are susceptible 
to couple with graphene states which makes the graphene monolayer on SiC (3 x 3) a nearly 
ideal system for investigating low energy excitations. 

Acknowledgments 

This work was supported by the French ANR ("GraphSiC" project N° ANR-07-BLAN- 
0161) and by the Region Rhone-Alpes ("Cible07" and "CibleOS" programs). F. Hiebel held 
a doctoral fellowship from la Region Rhone-Alpes. 



^ A. K. Geim and K. S. Novoselov, Nature Mater. 6, 183 (2007). 

2 A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, Rev. Mod. 
Phys. 81, 109 (2009). 

^ K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, M. I. Katsnelson, I. V. Grigorieva, S. V. 

Dubonos, and A. A. Firsov, Nature 438, 197 (2005). 
* Y. Zhang, Y. W. Tan, H. L. Stormer, and P. Kim, Nature 438, 201 (2005). 

5 M. I. Katsnelson, K. S. Novoselov, and A. K. Geim, Nat. Phys. 2, 620 (2006). 

6 A. F. Young and P. Kim, Nat. Phys. 5, 222 (2009). 

X. Wu, X. Li, Z. Song, C. Berger, and W. A. de Heer, Phys. Rev. Lett. 98, 136801 (2007). 

15 



^ E. McCann, K. Kechedzhi, V. I. Fal'ko, H. Suzuura, T. Ando, and B. L. Altshuler, Phys. Rev. 
Lett. 97, 146805 (2006). 

^ K. Bolotin, K. Sikes, Z. Jiang, M. Klima, G. Fudenberg, J. Hone, P. Kim, and H. Stormer, Solid 

State Commun. 146, 351 (2008). 

X. Du, I. Skachko, A. Barker, and E. Y. Andrei, Nature Nanotechnology 3, 491 (2008). 
1^ K. I. Bolotin, K. J. Sikes, J. Hone, H. L. Stormer, and P. Kim, Phys. Rev. Lett. 101, 096802 
(2008). 

12 V. V. Cheianov, V. Fal'ko, and B. L. Altshuler, Science 315, 1252 (2007). 

1^ L. A. Ponomarenko, F. Schcdin, M. I. Katsnelson, R. Yang, E. W. Hill, K. S. Novoselov, and 
A. K. Geim, Science 320, 356 (2008). 

1^ K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. Zhang, S. V. Dubonos, I. V. Grig- 
orieva, and A. A. Firsov, Science 306, 666 (2004). 

1^ I. Meric, M. Y. Han, A. F. Young, B. Ozyilmaz, P. Kim, and K. L. Shepard, Nature Nanotech- 
nology 3, 654 (2008). 

16 F. Schedin, A. K. Geim, S. V. Morozov, E. W. Hill, P. Blake, M. I. Katsnelson, and K. S. 

Novoselov, Nature Mater. 6, 652 (2007). 
1^ A. Griineis and D. V. Vyalikh, Phys. Rev. B 77, 193401 (2008). 

1^ T. Brugger, S. Giinthcr, B. Wang, H. Dil, M.-L. Bocquet, J. Osterwalder, J. Wintterlin, and 

T. Greber, Phys. Rev. B 79, 045407 (2009). 
19 P. Sutter, M. S. Hybertsen, J. T. Sadowski, and E. Sutter, Nano Lett. 9, 2654 (2009). 
2° C. Berger, Z. Song, T. Li, X. Li, A. Y. Ogbazghi, R. Feng, Z. Dai, A. N. Marchenkov, E. H. 

Conrad, P. N. First, et al., J. Phys. Chem. B 108, 19912 (2004). 

21 C. Berger, Z. Song, X. Li, X. Wu, N. Brown, C. Naud, D. Mayou, T. Li, J. Hass, A. N. 
Marchenkov, et al.. Science 312, 1101 (2006). 

22 W. A. de Heer, C. Berger, X. Wu, P. N. First, E. H. Conrad, X. Li, T. Li, M. Sprinkle, J. Hass, 
M. L. Sadowski, et al.. Solid State Commun. 143, 92 (2007). 

2^ A. J. van Bommel, J. E. Crombeen, and A. van Tooren, Surf. Sci. 48, 463 (1975). 

24 I. Forbeaux, J.-M. Themlin, and J.-M. Debever, Phys. Rev. B 58, 16396 (1998). 

25 L Forbeaux, J. M. Themlin, and J. M. Debever, Surface Science 442 (1999). 

26 S. Kim, J. Ihm, H. J. Choi, and Y.-W. Son, Phys. Rev. Lett. 100, 176802 (2008). 

2^ F. Varchon, P. Mallet, J.-Y. Veuillen, and L. Magaud, Phys. Rev. B 77, 235412 (2008). 

16 



28 A. Mattausch and O. Pankratov, Phys. Rev. Lett. 99, 076802 (2007). 

29 F. Varchon, R. Feng, J. Hass, X. Li, B. N. Nguyen, C. Naud, P. Mallet, J.-Y. Veuillen, C. Berger, 
E. H. Conrad, et al., Phys. Rev. Lett. 99, 126805 (2007). 

K. V. Emtscv, F. Speck, T. Seyller, L. Ley, and J. D. Riley Phys. Rev. B 77, 155303 (2008). 

31 P. Mallet, F. Varchon, C. Naud, L. Magaud, C. Berger, and J.-Y. Veuillen, Phys. Rev. B 76, 
041403(R) (2007). 

32 C. Riedl, U. Starke, J. Bernhardt, M. Franke, and K. Heinz, Phys. Rev. B 76, 245406 (2007). 

33 G. M. Rutter, N. P. Guisinger, J. N. Grain, E. A. A. Jarvis, M. D. Stiles, T. Li, P. N. First, 
and J. A. Stroscio, Phys. Rev. B 76, 235416 (2007). 

3^ P. Lauffer, K. V. Emtsev, R. Graupner, T. Seyller, L. Ley, S. A. Reshanov, and H. B. Weber, 
Phys. Rev. B 77, 155426 (2008). 

35 A. Bostwick, T. Ohta, T. Seyller, K. Horn, and E. Rotenberg, Nat. Phys. 3, 36 (2007). 

36 S. Y. Zhou, G.-H. Gweon, A. V. Fedorov, P. N. First, and W. A. de Heer et al.. Nature Mater. 

6, 770 (2007). 

3^" E. Rotenberg, A. Bostwick, T. Ohta, J. L. McChesney, T. Seyller, and K. Horn, Nature Mater. 

7, 258 (2008). 

38 S. Zhou, D. Siegel, A. Fedorov, F. Gabaly, A. Schmid, A. C. Neto, D.-H. Lee, and A. Lanzara, 
Nature Mater. 7, 259 (2008). 

39 V. W. Brar, Y. Zhang, Y. Yayon, T. Ohta, J. L. McChesney, A. Bostwick, E. Rotenberg, 
K. Horn, and M. F. Crommie, Appl. Phys. Lett. 91, 122102 (2007). 

^° I. Brihuega, P. Mallet, C. Bena, S. Bose, C. Michaelis, L. VitaH, F. Varchon, L. Magaud, 
K. Kern, and J. Y. Veuillen, Phys. Rev. Lett. 101, 206802 (2008). 

41 F. Hiebel, P. Mallet, F. Varchon, L. Magaud, and J.-Y. Veuillen, Phys. Rev. B 78, 153412 

(2008) . 

42 U. Starke and C. Riedl, J. Phys.: Condens. Matter 21, 134016 (2009). 

43 L. Magaud, F. Hiebel, F. Varchon, P. Mallet, and J.-Y. Veuillen, Phys. Rev. B 79, 161405(R) 

(2009) . 

44 A. L. Vazquez de Parga, F. Calleja, B. Borca, J. M. C. G. Passeggi, J. J. Hinarejos, F. Guinea, 
and R. Miranda, Phys. Rev. Lett. 100, 056807 (2008). 

45 J. Hass, R. Feng, J. E. Millan-Otoya, X. Li, M. Sprinkle, P. N. First, W. A. de Heer, E. H. 
Conrad, and C. Berger, Phys. Rev. B 75, 214109 (2007). 

17 



A. Varykhalov, J. Sanchez-Barriga, A. M. Shikin, C. Biswas, E. Vescovo, A. Rybkin, 
D. Marchenko, and O. Rader, Phys. Rev. Lett. 101, 157601 (2008). 
4^ H. E. Hoster, M. A. Kulakov, and B. BuUemer, Surf. Sci. Lett. 382, 658 (1997). 

I. Amidror, The Theory of the Moire Phenomenon (Kluwcr Academic Publishers, 1999). 

49 W.-T. Pong and C. Durkan, J. Phys. D: Appl. Phys. 38, R329 (2005). 

50 F. Varchon, P. Mallet, L. Magaud, and J.-Y. Veuillen, Phys. Rev. B 77, 165415 (2008). 

5^ A. T. N'Diaye, J. Coraux, T. N. Plasa, C. Busse, and T. Michely, New J. Phys. 10, 043033 
(2008). 

52 S. Marchini, S. Giinthcr, and J. Wintterhn, Phys. Rev. B 76, 075429 (2007). 

53 M.-H. Tsai, C. S. Chang, J. D. Dow, and L S. T. Tsong, Phys. Rev. B 45, 1327 (1992). 

54 Occasionally, a small contribution from the previously quoted in-gap states of the SiC(3 x 3) is 
seen in the range [— l.OF, — 0.5F]. 

55 We note that the measured decrease rate of the corrugation outside the gap depends on the 
tip state. These variations can be attributed to a varying sensitivity of the tips to interface or 
graphene states. A weak graphene signal superimposed on the substrate-related structures is 
sometimes observed up to 0.5 eV beyond the gap. 

56 S. Pons, P. Mallet, L. Magaud, and J. Y. Veuillen, Europhys. Lett. 61, 375 (2003). 

57 E. Sutter, D. P. Acharya, J. T. Sadowski, and P. Sutter, Appl. Phys. Lett. 94, 133101 (2009). 
5^ B. Wang, M.-L. Bocquet, S. Marchini, S. Giinther, and J. Wintterlin, Phys. Chem. Chem. Phys. 

10, 3530 (2008). 

59 D. Martoccia, P. R. Willmott, T. Brugger, M. Bjorck, S. Giinther, C. M. Schlepiitz, A. Cervel- 
lino, S. A. PauH, B. D. Patterson, S. Marchini, et al., Phys. Rev. Lett. 101, 126102 (2008). 

6° A. B. Preobrajenski, M. L. Ng, A. S. Vinogradov, and N. Martensson, Phys. Rev. B 78, 073401 
(2008). 

6^ B. Wang, M.-L. Bocquet, S. Giinther, and J. Wintterlin, Phys. Rev. Lett. 101, 099703 (2008). 

62 L. Vitah, C. Riedl, R. Ohmann, L Brihuega, U. Starke, and K. Kern, Surf. Sci. 602, L127 
(2008). 

63 C.-H. Park, L. Yang, Y.-W. Son, M. L. Cohen, and S. G. Louie, Nature Mater. 4, 213 (2008). 

64 L Pletikosic, M. Kralj, P. Pervan, R. Brako, J. Coraux, A. T. N'Diaye, C. Busse, and T. Michely, 
Phys. Rev. Lett. 102, 056808 (2009). 



18 



d) « 

5 

£ ' 

c 

— 3 
0. 

2 
1 



5 10 15 20 25 30 

a (deg) 

FIG. 1: (Color online) (a) 50 x 50 nm^ STM image of G_(3 x 3) islands on SiC(OOOT) with 
enhanced contrast to reveal their super lattices. Sample bias: —2.5 V. (b) LEED pattern of the 
sample, primary energy: 78 eV. Circles indicate first order SiC(l x 1) and SiC(3 x 3) spots. The 
relevant high order SiC(3 x 3) spots close to the graphitic ring-shaped signal are indicated by 
arrows and their coordinates are given, (c) Schematic picture of the LEED pattern in (b). The 
dashed curve stands for the graphitic signal. The dashed arrow is the reciprocal lattice vector 
that corresponds to a graphene layer rotated by an angle a with respect to the SiC(l x 1) surface. 
The dotted one corresponds to the (2,2) reciprocal vector of the SiC(3 x 3) reconstruction. The 
reciprocal lattice vector of the moire pattern k]vi(a) constructed on these vectors is represented by 
a solid vector, (d) The calculated moire periods P{a) for (4,0), (3,1) and (2,2) SiC(3 x 3) reciprocal 
lattice vectors are represented in dotted, solid and dashed line respectively. Triangles correspond 
to measurements on STM images. 
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FIG. 2: (Color online) (a), (b): 8x8 nm^ STM images of G_3 x 3 islands with atomic resolution 
on the graphene layer. Two different moire patterns are observed: (a) shows a superstructure with 
"ball-like" contrast (P= 2.9 nm), typical of a moire pattern constructed on the (2,2) SiC(3 x 3) 
Fourier component, (b) shows a superstructure with "hole-like" contrast (P= 3.5 nm), typical of 
a moire pattern constructed on the (3, 1) SiC(3 x 3) Fourier component, (c), (d): 5x5 nm^ STM 
images of two other G_3 x 3 islands showing the same type of moire pattern as (a) and (b) resp. 
((c) P= 2.9 nm (d) P= 4.7 nm). The sample bias is chosen so that both graphene and SiC(3 x 3) 
lattices are detected and their stacking is thus visible, (e), (f): Schematic representations of the 
stacking of the graphene atoms on the substrate reconstruction deduced from images (c) and (d) 
respectively (the link between the STM images and the schematics is indicated by arrows). On 
(e), the Wigner-Seitz (pseudo) common cell is represented in solid line. On both illustrations, the 
dotted circle separates regions of the moire cell with different types of stacking. 
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FIG. 3: (Color online) 8x8 nm^ STM images of (a) an a = 30° G_3 x 3 island; sample bias: 
—1.5 V. The SiC(6\/3 x 6\/3)i?30° common cell is represented with dashed lines. Insert: 3 x 3nm^ 
STM image of the same island; sample bias: +100 mV.The graphene layer shows an AB symmetric 
honeycomb contrast, (b) an a = 29° G_3x3 island; sample bias: -1.65 V. The SiC(6\/3x6\/3)fl30° 
cell represented by dashed lines clearly does not coincide with the moire pseudo unit cell represented 
by solid lines. 
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FIG. 4: (Color online) (a) Inset: 20 x 20 nm^ STM image of a region with the bare SiC(3 x 3) 
reconstruction (right side) and a G_3 x 3 island (left side). Sample bias: —2.5 V. The I(V) curves 
are spatially averaged over the boxed regions (300 points) on the bare SiC(3 x 3) reconstruction 
(gray line) and on the G_3 x 3 island (black line) and over 15 points of the island edge (dashed line). 
For STS measurements: Setpoint voltage: —2.5 V; Setpoint current: 1.0 nA. (b) Corresponding 
dl/dV curves for the G_3 x 3 island (black line) and the bare SiC(3 x 3) reconstruction (gray line). 
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FIG. 5: (Color online) (a) A series of 8 x 6 nni'^ STM images taken at the same spot of a G_3 x 3 
island (a = 16°, P = 4.5 nm) at various sample biases (V^) (Common height range AZ = lA). The 
relative amplitude of the signals arising from the SiC(3 x 3) reconstruction and from the graphene 
layer changes with sample bias, (b) Laterally averaged profile taken over the boxed region in (a) 
for Vs = +10mV. (c) The peak-to-peak moire corrugation amplitude as a function of the sample 
bias voltage, on two similar G_3 x 3 islands (a = 16°, P = 4.5nm) of two different samples obtained 
with different tips. The island referred to as "island 1" is the island presented in (a). A dl/dV 
curve acquired on the bare SiC(3 x 3) reconstruction is also represented on the graph (solid line) 
in order to locate the substrate surface bandgap. The moire corrugation and the LDOS of the 
substrate show complementary behaviors. 
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FIG. 6: (Color online) Inset: 10 x 10 nm^ STM image of a G_3 x 3 island presenting a "ball- 
like" moire contrast, (a = 28°, P = 2.9 nm). Sample bias: —1.0 V. dl/dV curves are spatially 
averaged on 24 points of highest topographic regions (black) and of lowest topographic regions 
(gray). Setpoint voltage: —1.0 V; Setpoint current: 1.0 nA. No spatial variation in the electronic 
structure of graphene is observed in the voltage range [—1.0 V, +0.15 V]. 
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